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Abstract 

Synchronization of forced reactively coupled van der Pol oscilla¬ 
tors is investigated in the phase approximation. We discuss essen¬ 
tial features of the reactive coupling. Bifurcation mechanisms for the 
destruction of complete synchronization and possible quasi-periodic 
regimes of different types are revealed. Regimes when autonomous 
oscillators demonstrate frequency locking and beating regimes with 
incommensurate frequencies are considered and compared. 
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1 Introduction 

Synchronization is a common phenomenon in physics, chemistry, biology and 
other fields of science and technology DEI- Nowadays, different examples 
and variety of synchronization in ensembles of oscillators are topics of the 
research interest [3j. It has been found that systems containing even a small 
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number of elements demonstrate a great variety of different types of dy¬ 
namics m- However, a complete picture of synchronization is still not 
constructed. One of the fundamental problems is a problem of influence of 
an external signal to coupled oscillators. Some new aspects of this problem 
have been recently revealed for dissipative coupling. A corresponding phase 
equation which generalizes classical Adler equation has been obtained in [Dj. 
A case when two autonomous oscillators are locked has been investigated. 
A complete synchronization region for two oscillators with an external force 
has been detected, and degenerated saddle-node bifurcation of equilibrium 
states responsible for the destruction of such synchronization has been de¬ 
scribed. It has been shown that three-frequency quasi-periodicity appears 
in the phase model through a saddle-node bifurcation of invariant curves. 
External force synchronization has been investigated in |7j in situation when 
a locking regime in system of autonomous oscillators changes to a beating 
regime with incommensurable frequencies. An influence of asymmetry of 
coupling to synchronization picture has been investigated in [8j. Excitation 
of a chain of three phase oscillators with dissipative coupling has been con¬ 
sidered in [9]. Synchronization mechanisms beyond the phase approximation 
have been discussed in [ID]. A possibility of both the regimes observed in the 
phase model and a new kind of synchronization via suppression of oscillations 
has been shown. 

These results, however, refer to the simplest case when the coupling is 
dissipative. There is another type of coupling, namely reactive (or conser¬ 
vative) coupling. This is a situation when the coupling is realized directly 
through values of the variables and not through their rates mini. Reactive 
coupling occurs in many interesting problems. An actual example of such 
coupling is ion traps [12] . In such traps ions are confined using variable mi¬ 
crowave fields, which restrict a magnitude of radial oscillations of the ions 
and a constant electric field limiting the axial motions. Laser radiation of 
different frequencies provides instability and dissipation into the system. Re¬ 
active type of the coupling is caused by Coulomb repulsion between the ions. 
Moreover a reactive coupling arises in the array of coupled nanomechanical 
cavities as well as in some optical problems [ ID] , 

From a theoretical point of view, the reactive coupling is an essentially 
more delicate phenomenon than the dissipative coupling mini. So in the 
simplest case of two coupled oscillators, equations in the phase approximation 
can be obtained only if we deal with second order of smallness in the coupling 
parameter. For this, in contrast to the dissipative coupling, it is necessary to 


2 


take into account deviation of their orbits from the unperturbed states Ell- 
Phase bistability is possible for the reactive coupling when synchronization 
between oscillators may be in-phase as well as anti-phase. At the same time 
a problem of forced synchronization of reactively coupled oscillators has not 
yet been fully investigated. We should note the work [T3j having a style of the 
computer experiment. This work has demonstrated a possibility of two- and 
three-frequency quasi-periodicity and chaos in such a system. In the more 
recent paper [15], a structure of the parameter plane ’’frequency - amplitude 
of an external signal” has been investigated, but the most significant and 
complex region of small signal amplitudes has not been considered. At the 
same time, a phase model has not been constructed and the authors have 
investigated only the case when autonomous oscillators are locked. These 
works describe the variety of possible effects and bifurcations in such system 
not well enough. Therefore, we discuss in this paper the following questions: 

• How to write the correct phase equations for the forced system in case 
of the reactive coupling between oscillators? 

• What is the structure of complete synchronization region for such forced 
oscillators, and what bifurcation is responsible for the destruction of 
this region? 

• What effects do occur in the presence of phase bistability in case of an 
external forcing? 

• How are the various two-frequency and three-frequency regimes embed¬ 
ded in the parameter space? 

An important moment is a parallel consideration of two possible situ¬ 
ations. The first one is when autonomous oscillators are locked, and the 
second one is when they demonstrate oscillations in the beating regime. 
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2 Phase equations of reactively coupled os¬ 
cillators 


Consider a system of two reactively coupled van der Pol oscillators excited 
by an external harmonic signal: 

x — (A — x 2 )x + (1 ——)x + e{x — y) = B sin ojt , 

( 1 ) 

V ~ (A - y 2 )y + (1 + —)y + e(y - x) = 0. 

Here A is the control parameter, A is the frequency detuning between os¬ 
cillators, e is the parameter of reactive coupling, B is the signal amplitude, 
and oj is the signal frequency. Central frequency of the oscillators is taken as 
unity. 

We assume that control parameter, coupling value and signal frequency 
detuning from the unit frequency are small. In this case we can use the 
method of slowly varying amplitudes [I]. We represent the dynamic variables 
as 

x = a(t)e luJt + a*(t)e~ lult , y = b(t)e lut + b* . (2) 

and apply standard additional conditions: 

a(t)e iut + a*(t)e~ iut = 0, b{t)e iujt + b*{t)e~ iujt = 0. (3) 


where a(t) and b(t) are the complex amplitudes. Let us substitute the re¬ 
lations (j2J), (|3]) in the equations (0Q), multiply the obtained expressions by 
e~ lult and average the result for eliminating rapidly oscillating terms. Also 

B(e iuit - Be~ iut ) 

we use the relation Bsinut = -. After some mathematical 
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transformations we obtain the Landau-Stuart equations: 


1 2 A 
I.--- or- 

O cy 

2a = Xa + i - —a — \a\ a + ie(a — b) — 

OJ 


B 

2u? 


1 - oj 1 + 


A 


2b = Xb — i- 


OJ 


-b — \b\ 2 b + ie{b — a). 


( 4 ) 
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Now we introduce the frequency detuning of the signal from the central 

1 2 A 

A 

-— « —2fl=F—. Fur- 

u 2 


frequency: uj — 1 + fh Since u ~ 1, we can set 


thermore, the parameter A in Eqs.(J5J) can be eliminated by renormalization: 

t 12 A e B k , . 

t —> —, iz —> —, A —y —, e —y — and £>—>—. As a result, we obtain: 

A A A A A 


2 d = a + z (20 + — )a — \a\ 2 a + ie(a — b) ——, 
2 b = b — i(2Q — —)b — \b\ 2 b + ie{b — a). 


(5) 


Let us introduce the real amplitudes ri )2 and phases of the oscillators with 
respect to the external signal - 0 li2 using the relations a = Tqe l ^ 1 and b = 
r 2 e 1 ^ 2 . Then Eqs.(E]) yield: 


2 ri = ?’i — rj — er 2 sin 6 -cos i/q, 


2 r 2 = r 2 — rg + eri sin 0 , 

A r 2 5 

2ipi = —20 -he- 1 e cos 0 H-sin ipi , 

2 ri 2 

2 i /^ 2 = —20 + A + e — — e cos 
2 r 2 


( 6 ) 


Here 9 = ifti — ^ 2 is the relative phase of the oscillators. 

Let us now to take into account the deviation of orbits from stationary 
unperturbed values r x = r 2 = 1. For this we set r 1)2 = 1 + fi )2 , where f 1)2 are 
the small perturbations, and rewrite first and second equations in Eqs. (TB1) as: 


2ri = — 2f\ — £ sin 9 ——cos i/q, 2r 2 = —2r 2 + esind. (7) 

In these equations the amplitude perturbations are strongly damped, so there 
is rapid convergence to the stationary states: 


f\ — — -sin#—— cos " 01 , f 2 = ^ sin 9. 


( 8 ) 
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Substituting these expressions in third and fourth equations in Eqs.flH]) 
and retaining only the first and second order terms of the coupling parameter 
e gives: 


0i = — —— 0 + ^(1 — cos#) —— sin 26 + — sin-0!, 

(9) 

A £ £ 2 

02 =-H —(1 — cos 6) H-sin 26. 

4 2 4 

Here we assume that the amplitude B is of order e 2 and neglect the term Be. 
Physical explanation is that the external signal perturbs only phases of the 
oscillators and practically does not perturb amplitudes. Subtracting Eqs.(j9]) 
from each other, we obtain: 


6 = - 

01 = 


A 
~2 ~ 

A 

-J 


— sin 26 H-sin0i, 

2 4 ^ ’ 

£ £ 2 B 

-OH— (1 — cos 6) -sin 2# + — sin 0!. 

2 4 4 


( 10 ) 


Eqs. (lT0D are the desired phase equations for forced reactively coupled 
oscillators. In the absence of coupling £ = 0 we obtain phase equation of 
excitation of the first oscillator [1], On the other hand, in the absence of 
external signal B = 0 Eqs. fllO]) can be written as the well-known equation 
for the relative phase 6 = 0i — 0 2 of two reactively coupled oscillators HU: 
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A 

~2 


— sin 2 6. 
2 


( 11 ) 


This equation describes an ability of the beating regime at |A| > e 2 and 
mutual locking of oscillators at |A| < e 2 . We can see also from Eq. fTTT]) that 
two oscillators may exhibit both the in-phase (6 ~ 0) and anti-phase (6 fa 7r) 
stable synchronization. There are also two unstable synchronous regimes in 
the intermediate values of 6 HU- 

Model (HOD was obtained by a number of restrictions on the parameters. 
However, it turns out, as we will show, quite effective in a wider area. 


6 


3 Analysis of equilibrium states. Saddle-node 
bifurcations 


One of the basic regimes in the forced system is a complete synchronization 
when an external signal is locking for both oscillators. In this case relative 
phases of the oscillators with respect to the external signal ipi and ij > 2 are 
constant in time. This regime corresponds to a stable equilibrium state of 
the system (fTOlh Let us find the region of complete synchronization on the 
parameter plane and discuss the mechanisms of its destruction. 

In the case of dissipative coupling, one stable node, one unstable node 
and two saddles coexist in the phase model [6]. Destruction of the complete 
synchronization is caused by the specific variant of saddle-node bifurcation 
when all the four equilibrium states merge in pairs simultaneously and dis¬ 
appear, and stable and unstable invariant curves occur from their manifolds 
[6]. In the case of reactive coupling, the mechanisms of synchronization and 
its destruction are changing and become more complicated. 

To find the equilibrium states, we set i\)\ = ^2 = 0 in Eqs. ffTOlk 

A 1 9 

—-— -e sin 29 + b sin '0i = 0, 

( 12 ) 

—D — — + ^(1 — cos 9) — -e 2 sin 29 + bsm'ipi = 0. 


Here b 


B 

T 


The perturbed matrix of the system (TTOl) is 


M 


( —£ 2 cos 29 

£ £ 2 

- sin 9 -cos 29 

V 2 2 


b cos V’i 
b cos ijj 1 


(13) 


The condition of saddle-node bifurcation is zero Jacobian of this matrix, i.e. 
det M = cos ifti [e cos 29 + sin 9] = 0. The obtained equation can be split into 
two equations: 

e cos 20 + sind = 0, cos-i/q = 0. (14) 

The Erst one leads to the relation 


sin 9 


1 =F \/l + 8e 2 
4e 


(15) 
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At £ < 1 it has a solution only if we use an upper sign. Subtracting one equa¬ 
tion from another in Eqs. fll2D gives the following conditions for the saddle- 
node bifurcation: 


/\ £ 2 . 

Vt =-b -(1 — cos#) + -£ 2 sin2#, 

4 2 4 


„ • ,1-VT+8? N 

#i = arcsm(-), 

4s 


6o = it — arcsim 


1 - Vl + 8e 2 , 

' 4s ' 


(16) 


Conditions (fTTfl) dehne the locking band bounded by two vertical lines on 
the (Jl, b ) parameter plane, so that Q mm < hi < Q max (see Fig.la). Width 
of the locking band is determined by the coupling parameter. It is easy to 
show that at small coupling values it has a width approximately equal to e. 



Figure 1: Lines (fTTfl) and (fTTl) for e = 0.3; a) beating regime of autonomous 
oscillators, A = 0.3; b) locking regime of autonomous oscillators, A = 0.05. 

and are the eigenfrequencies of the in-phase and out-phase modes. 
Unsigned arrows indicate locking frequencies in the unstable regimes. 

Now consider the second condition cos'i/q = 0. Then sin yjj = ±1 and 
Eqs.(H2D, (HI give 

2b = ±(A + s 2 sin 26), U = — + ^(1 — cos 9) + -s 2 sin 20. (17) 

The lines have a different configuration depending on whether the 
autonomous oscillators are in beating (A > s 2 ) or in locking (A < 2 ) regimes. 
















Both of these cases are shown in Fig.l on the left and right side, respectively. 
Here and below, we need to select the normalized values of the coupling 
parameter e. We use the value e = 0.3. The smallness of the parameter e 2 
(in this case e 2 = 0.09) is a condition of the effectiveness of the phase model. 

At the same time, the choice of e = 0.3 will further demonstrate certain 
differences of the original system from the phase model. 

In the case of beating regime A > e 2 the condition b ^ 0 is always 
satisfied. Therefore, the line (1171) does not touch frequency axis and the 
complete synchronization region has an amplitude threshold in the beating 
regime of autonomous oscillators. This line forms a closed curve in the form 
of the ’’eight” number, Fig.la. It is easy to show that the vertical lines (fTUD 
touch this curve. 

In the case of locking regime of autonomous oscillators the configuration of 
lines becomes more complicated, Fig.lb. Now they have four characteristic 
tips at b = 0 marked by arrows in Fig.lb. If the signal amplitude tends 
to zero b —> 0, Eqs. flTTD give the conditions for locking frequencies of two 
autonomous reactively coupled oscillators. The values marked in Fig. lb by 
arrows correspond to the eigenfrequencies of these modes. VL* n and Q* out 
indicate the frequencies of stable in-phase and out-phase regimes, and similar 
points for unstable regimes are marked by unsigned arrows. Thus, the in- 
phase and out-phase mode arc locked at different frequencies. 

4 Complete picture of bifurcations. Andronov- 
Hopf bifurcation 

However, it turns out that there are two different types of locking of both 
oscillators by an external signal in a system with reactive coupling. The 
Erst one is defined by a stable equilibrium of the system (TTOh and has been 
discussed above. The corresponding phase portrait is presented in Fig.2, 
fragment 1. The second type of locking corresponds to a stable limit cycle 
of the system (HUD shown in Fig.2, fragment 6. In this locking regime, the 
phases of oscillators are not constant and oscillate near some equilibrium 
value. Therefore, even in the phase approximation, a new line defined by the 
period of phase oscillations appears in the spectrum of oscillators. 

From this it follows that the Andronov-Hopf bifurcation (the stable limit 
cycle birth bifurcation) is possible in the phase system with the reactive 


9 


coupling. Therefore, we consider a more complete bifurcation picture of the 
analyzed system shown in Fig. 3 and revealed numerically]]] This figure indi¬ 
cates also two cases: beating and locking regimes of autonomous oscillators. 
For the beating regime (Fig.3a), there are two branches of saddle-node bi¬ 
furcation lines SN revealed in the analytical study. They touch each other 
at degenerate cusp points DC. The saddle-node bifurcation lines contain 
the segments SNi where merging of stable equilibrium and saddle occurs 
(solid line) and the segments SN 2 where merging of unstable equilibrium 
and saddle takes place (dashed line). The Andronov-Hopf bifurcation H is 
also possible. Supercritical bifurcation Hi is responsible for emergence of a 
stable limit cycle, and subcritical bifurcation H 2 for is responsible disappear¬ 
ance of an unstable limit cycle. The Andronov-Hopf bifurcation lines have 
common points with the saddle-node bifurcations line, the Bogdanov-Takens 
points BT. Grey color in Fig.3 shows the regions where the system has at 
least one equilibrium state. 

Examples of the phase portraits at characteristic points labeled with num¬ 
bers in Fig.3 are given in Fig.2. The simplest mechanism of the complete 
synchronization destruction is concerned with the saddle-node bifurcation of 
equilibrium states when the stable equilibrium and the saddle merge and 
an invariant curve responsible for the quasi-periodic oscillations appears (see 
the transition from point 1 to point 4 in Figs.2,3). Unlike the case of dissi¬ 
pative coupling [5], the second pair of equilibrium points continues to exist. 
Decreasing of the signal amplitude leads to their mergence but on the lower 
branch of the saddle-node bifurcation line SN 2 . 

The second mechanism becomes apparent if we move from point 2 to 
point 4 through point 3. It consists of the following. Initially, an unstable 
limit cycle appears as a result of the nonlocal bifurcation. This limit cycle 
converges to a stable equilibrium point which disappears through the inverse 
Andronov-Hopf bifurcation H 2 . 

Let us consider a more complex situation when the locking of autonomous 
oscillators is observed (Fig.3b). I 11 this situation, another mechanism for the 
destruction of synchronization is possible (see the transition from point 5 
to points 6 and 7). Amplitude decreasing of an external signal leads now 
to the supercritical Andronov-Hopf bifurcation on the line H\. As a result, 
stable focus becomes unstable and gives birth to stable limit cycle in the 

^^Note that analytical investigation in Sec.2 is an essential supplement to numerical 
analysis due to the complexity of the bifurcation picture. 
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Figure 2: Phase portraits of the system (II Oh . Numbers correspond to the 
points in Fig.3. Values of the parameters are A = 0.3, 12 = 0.21 and B = 0.8 
(Fragment 1); A = 0.3, 12 = 0.21 and B = 0.75 (Fragment 2); A = 0.3, 
12 = 0.22 and B = 0.67 (Fragment 3); A = 0.3, 12 = 0.17 and B = -0.66 
(Fragment 4); A = 0.05, 12 = 0.15 and B = 0.5 (Fragment 5); A = 0.05, 
12 = 0.15 and B = 0.355 (Fragment 6); A = 0.05, 12 = 0.16 and B = 0.23 
(Fragment 7); A = 0.05, 12 = 0.29 and B = 0.7 (Fragment 8); A = 0.05, 
12 = 0.31 and B = 0.035 (Fragment 9). 


system. Furthermore, unstable limit cycle is born from a separatrix loop. (In 
fragment 6 this cycle is almost touching the separatrix.) Stable and unstable 
limit cycles come closer to each other and then merge and disappear on the 
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Figure 3: Bifurcation lines and points of the system (HU on the (hi, B) 
parameter plane for e = 0.3; a) beating regime, A = 0.3; b) locking of au¬ 
tonomous oscillators, A = 0.05. SN is the saddle-node bifurcation of equi¬ 
librium states, Hi is the supercritical Andronov-Hopf bifurcation, Ho is the 
subcritical Andronov-Hopf bifurcation, SNC is the saddle-node bifurcation 
of limit cycles, DC is the degenerate cusp point, BT is the Bogdanov-Takens 
point, GH is the generalized point of the Andronov-Hopf bifurcation. 


saddle-node bifurcation line of limit cycles SNC. Remaining in the phase 
space unstable focus is transformed into an unstable node, merges with saddle 
and disappears on the line SN. 

Note that for the two-parameter analysis the number of variants is suffi¬ 
ciently high and depends on the path in the parameter plane. At the same 
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time, possible scenarios depend on the location of codimension-two points, 
the Bogdanov-Takens points BT, and generalized points of the Andronov- 
Hopf bifurcation GH , in which a supercritical Andronov-Hopf bifurcation 
turns into a subcritical one. The saddle-node bifurcation line of limit cycles 
SNC comes also into the point GH. 

In the considered situation, there are also new variants of the phase mul¬ 
tistability. In the point 8, the stable invariant curve coexists with the stable 
equilibrium. Equilibrium state corresponds to the phase shift between oscil¬ 
lators equal approximately to it. In this case, out-phase oscillations of the 
oscillators are locked by an external signal. Stable quasi-periodic oscillations 
occur when the phases are approximately identical. Beyond the synchroniza¬ 
tion tongue (point 9 in Fig.2,3), there is a situation of coexistence of two 
stable invariant curves corresponding to the in-phase and out-phase locking 
of relative phase of oscillators. In similar way, in-phase oscillations of the 
oscillators are locked by an external signal near the first tip in Fig.3b. 

Therefore, the bifurcation picture of a system with the reactive coupling 
differs significantly from that of a system with the dissipative coupling and 
is more complicated. 


5 Different types of regimes in the parameter 
plane 

Now we consider a more complete picture of different regimes possible in the 
system co. For this, we supplement the results of bifurcation analysis by 
the method of the charts of Lyapunov exponents. This method is highly 
effective for analysis of the systems with multi-frequency dynamics [5]. We 
calculate both Lyapunov exponents A[ and A 2 of the system at each grid 
point on the parameter plane. Then we color these points in accordance with 
the values of Lyapunov exponents, so the following regimes are visualized: 

1. P is a periodic regime with Ai < 0 and A 2 < 0 (red color); 

2. T 2 is a two-frequency quasi-periodic regime with Ai — 0 and A 2 < 0 
(yellow color); 

3. T 3 is a three-frequency quasi-periodic regime with A 3 = 0 and A 2 = 0 
(blue color). 
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Obtained charts are given in Fig.4. Fig.4a shows the situation when 
oscillators demonstrate beating regime in the absence of an external signal. 
In this case one can see the system of two-frequency synchronization tongues 
with the tips on the O -axis. These tongues are located inside the region of 
three-frequency quasi-periodicity T 3 . Region of the periodic regimes P has 
the threshold by amplitude of the signal. 



Figure 4: Charts of Lyapunov exponents for the phase system ffTOT) : e = 0.3; 
a) A = 0.3; b) A = 0.05. 

Another situation, when autonomous oscillators are locked in the absence 
of an external signal, is presented in Fig.4b. Here, the region of the three- 
frequency regimes T 3 has the threshold by amplitude of the signal. Inside this 
region, there is a set of narrow two-frequency resonance tongues of different 
type. The region of the periodic regimes P touches the fl-axis in two points. 
These points correspond to the resonances of an external signal with in- 
phase and out-phase modes, see Fig. lb. For small amplitudes of the external 
signal, only two-frequency quasi-periodicity is observed. Note that the chart 
of Lyapunov exponents reveals fine complex structure of the two tips in 
Fig.4b which is due to the fact that the region of complete synchronization 
is bounded by the saddle-node and Andronov-Hopf bifurcation lines. 
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6 Dynamics of the initial system 

Let us give some illustrations for the initial system (JT]). Its properties de¬ 
pend on the control parameter A. We set A = 0.1 which is much less than 
one. Therefore, one can expect a good agreement with the results for the 
phase model. On the other hand, this value is not very small, and possible 
differences may be observed. Corresponding chart of Lyapunov exponents 
is presented in Fig.5. Besides the regimes which are typical for the phase 
model, there are colored black chaotic regimes C with positive largest Lya¬ 
punov exponent. 



Figure 5: Charts of Lyapunov exponents for the system (Jl]); A = 0.1, £ = 0.03; 
a) A = 0.03; b) A = 0.005. If compare to Fig.4, all the parameters are 
normalized by A. 

Comparison of Fig.4 and Fig.5 allows to conclude that the phase model 
describes the main elements of synchronization picture sufficiently well even 
at A = 0.1. Fig.4a and Fig.5a demonstrate not only the same configuration of 
the region of complete synchronization P but also identical structure of the 
two-frequency tongues. At the same time, the tips of the region of complete 
synchronization looks similar in detail in Fig.4b and Fig.5b constructed for 
the case when autonomous oscillators are locked. Certain differences occur 
with increase of the signal amplitude B. (This is obvious for assumptions 
made for the phase model.) Thus, chaos occurs in the overlapping area of two- 
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frequency tori (Fig.5a). The region of three-frequency regimes are replaced 
significantly by the region of two-frequency regimes in (Fig.5b). Note that 
the phase model will work better than less control parameter A and coupling 
constant e. 

Fig.6 shows projections of the attractors on the variable planes for the 
first and second oscillators of the system ([I]). There is also (x, y ) plane 
which is suitable for visual type definition of the regime. Fig.6a represents 
the two-frequency torus, Fig.6b shows the three-frequency torus, and Fig.6c 
illustrates the chaotic regime. One can see that the regimes are sequentially 
complicated. For the invariant tori, attractor portraits look like the orbits 
similar to circles. This indicates the phase approximation efficiency. For 
the chaotic regime, orbit of the first oscillator is strongly disturbed, and 
trajectory attends neighborhood of the origin. Now the phase model is not 
useful. 

Note that the phase model will work better than smaller coupling constant 
e. We will reduce tripled the coupling constant that would confirm this. The 
corresponded chart of Lyapunov exponents is presented in Fig.7. From the 
comparison of Figs. 7 and 4b, we see that in this case the phase model 
becomes practically precision. 

7 Conclusion 

When describing synchronization of two reactively coupled oscillators by an 
external signal, we need to consider effects which are of the second order 
in coupling parameter. Locking of both oscillators by an external signal 
can be accurate or not. When the locking is accurate, there is a stable 
equilibrium in the phase equations. For the non-accurate locking, there is a 
limit cycle on the phase plane. An external signal of small amplitudes can 
lock both the in-phase and out-phase oscillations of the oscillators. Basic 
bifurcation mechanisms of the synchronization destruction are as follows: 
merging of the stable and unstable equilibriums; subcritical Andronov-Hopf 
bifurcation which results in disappearing of the stable equilibrium; saddle- 
node bifurcation of the stable and unstable limit cycles. The bifurcation 
picture includes also the codimension-two bifurcations: Bogdanov-Takens 
and generalized Andronov-Hopf bifurcation points. A feature of the system 
of forced reactively coupled oscillators is an ability of the phase bistability 
which consists in coexistence of different types of stable regimes (equilibrium 
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Figure 6: Examples of attractor projections for the system of three excited 
oscillators (JH) ; A = 0.1, A = 0.03, e = 0.03. a) Two-frequency torus, O = 
0.975, B = 0.04; b) three-frequency torus, fl = 0.95, B = 0.01, c) chaotic 
attractor, Q = 0.975, B = 0.0375. 

points, limit cycles, invariant curves). 
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